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^ ' Abstract 

' Classical random walks on weU-behaved graphs are rapidly mixing towards the uniform distribution. 

|. Moore and Russell showed that a continuous quantum walk on the hypercube is instantaneously uniform 

' mixing. We show that the continuous-time quantum walks on other well-behaved graphs do not exhibit 

this uniform mixing. We prove that the only graphs amongst balanced complete multipartite graphs 
that have the instantaneous uniform mixing property are the complete graphs on two, three and four 
vertices, and the cycle graph on four vertices. Our proof exploits the circulant structure of these graphs. 
lO ' Furthermore, we conjecture that most complete cycles and Cayley graphs lack this mixing property as 

. well. 

o _ 

5^ ■ 1 Introduction 

o ■ 

CN ' Two pervasive algorithmic ideas in quantum computation are Quantum Fourier Transform (QFT) and am- 

plitude amplification (see [2|). Most subsequent progress in quantum computing owed much to these two 
beautiful ideas. But there are many problems whose characteristics matches neither the QFT nor the ampli- 
tude amplification mold (e.g., the Graph Isomorphism problem). This begs for new additional tools to be 
^ ■ discovered. 

^ , A natural way to discover new quantum algorithmic ideas is to adapt a classical one to the quantum 

qh| model. An appealing well-studied classical idea in statistics and computer science is the method of random 

walks 001. Recently, the quantum analogue of classical random walks has been studied in a flurry of 
works lEl |6l IHl |9| 20). The works of Moore and Russell ||9l and Kempe iflOl showed faster bounds 
on instantaneous mixing and hitting times for discrete and continuous quantum walks on the hypercube 
■ (compared to the classical walk). 

The focus of this note is on the continuous-time quantum walk that was introduced by Farhi and Gut- 
mann |5J. In a subsequent paper, Childs et al. j6l gave a simple example where the classical and (continuous- 
time) quantum walks exhibit a different behavior in hitting time statistics. The goal of this note is to show 
further qualitative differences (by elementary means) between the two models in their mixing time behav- 
iors. A main result is that, on the natural class of complete and balanced complete multipartite graphs, 
only the complete graphs Kn on n = 2, 3, 4 vertices and the cycle graph of size four (i.e., ^2,2) have the 
instantaneous uniform mixing property, i.e., there exists a time when the probability distribution function of 
the quantum walk is exactly the uniform distribution. This is distinctly different from the classical walk and 
from the discrete quantum walk. The proofs of these exploit heavily the circulant structure of these graphs. 

A recent work by Childs et al. 1 1 1 gave an interesting and powerful algorithmic application of continuous 
quantum walks. 
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Figure 1 : The only balanced multipartite graphs with the continuous mixing property. From left to right: 

K2,K3 = C3, Ki, and ^2,2 = C4. 



2 Continuous-time quantum walks 

Continuous-time quantum walks was introduced by Farhi and Gutmann f51 (see also Our treatment, 

though, follow closely the analysis of Moore and Russell |9| which we review next. Let G = {V, E) be a 
simple, undirected, connected, d-regular n- vertex graph. Let A be the adjacency matrix of G that is given 
as 

r 1 ii{3,k)eE 

~ \ otherwise ^ ' 

We define the transition matrix H = (treated as the Hamiltonian of the quantum system). Let the initial 
amplitude wave function of the particle be |^o) = |0). The the amplitude wave function at time t, is given 
by Schrodinger's equation, as 

ih^JiPt) = H\^t). (2) 



dt' 



or (assuming from now on that h = 1) 



|V^t) = e-^^*|V'o). (3) 

This is similar to the model of a continuous-time Markov chain in the classical sense (see fTT] ). It is more 
natural to deal with the Laplacian of the graph, which is defined as L = A — D, where D is a. diagonal 
matrix with entries Djj = deg{vj). This is because we can view L as the generator matrix that describes an 
exponential distribution of waiting times at each vertex. But on d-regular graphs, = ^1, and since A and 
D commute, we get 

g-itL ^ ^-it{A~^I) ^ ^-it/d^-itA 

which introduces an irrelevant phase factor in the wave evolution. 
The probability that the particle is at vertex j at time t is given by 

PtU) = \{Mt)\^- (5) 

Since H is Hermitian, the matrix Ut = e~*^* is unitary. If {Xj, \ zj))j are the eigenvalue and eigenvector 
pairs of H, then (e"*^^*, \zj))j are the eigenvalue and eigenvector pairs of Ut- Because H is symmetric, 
there is an orthonormal set of eigenvectors, say {\zj) : j G [n]} (i.e., H is unitarily diagonalizable). So, if 
iV'o) = then 

\i^t) = ^a,e-'^^%). (6) 

j 

Hence, in order to analyze the behavior of the quantum walk, we follow its wave-like patterns using the 
eigenvalues and eigenvectors of the unitary evolution Ut- To observe its classical behavior, we collapse the 
wave vector into a probability vector using Equation |5] 
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In the following, we define the notion of (instantaneous) uniform mixing property of quantum walks 
(defined earlier in l9l fT0l ). 

Definition 1 A graph G = {V, E) has the instantaneous uniform mixing property if there exists t G M"*", 
such that {ipt) of the continuous-time quantum walk on G satisfies Pt{j) = 1/\V\, for all j G V. 

A random walk on G is called simple if the transition probability matrix is H = ^A, where A is 
the adjacency matrix of the d-regular graph G. The walk is called lazy if, at each step, the walk stays 
at the current vertex with probability ^ and moves according to H with probability ^. For continuous- 
time quantum walks, these two notions are equivalent modulo a time scaling. The argument is similar to 
Equation |4] by exploiting the commutativity of A and / that introduces an irrelevant phase factor in the 
amplitude expression. 



3 A simple example 

This section describes a single example that illustrates differences between classical (discrete and continu- 
ous) random walk and continuous-time quantum walk. The example concerns a random walk on K2. Let 
the Pauli X matrix and the Hadamard matrix H be given as 

The adjacency matrix of K2 is given by A = X. Since K2 is bipartite, a classical simple random walk 
starting at the first vertex will oscillate between the two states. So the simple random walk will never reach 
the uniform distribution. However, the lazy random walk, with transition matrix ^{I + A), will converge to 
the uniform distribution after a single step. 

Let us consider a continuous-time walk with generator matrix 

This is viewed as a two-state chain where a transition is made after a waiting time that is exponentially 
distributed with rates a and /3 from states 1 and 2, respectively (see |11|, page 156). Note that Q has 
eigenvalues 0, — (a + P) with the following respective eigenvectors 



a)' (-1 



(9) 



Let P{t) be the matrix with entries Pjk{t) = Pr(X(t) = k \ X{0) = j) which describes the evolution of 
the continuous-time walk after t steps. It is known that 

Ap(t) = QP{t) (10) 

which implies P{t) = e**^. Note that Q can be diagonalized using the following matrix B and its inverse: 
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Therefore, 



So, the uniform distribution is reached in the Umit if and only if a = /5 = 1. 

Finally, we consider the continuous-time quantum walk on K2. The wave amplitude vector at time t is 
given by 




Thus, collapsing this wave vector to its probability form, the expression 

is obtained which reaches the uniform distribution periodically at t = fevr + 7r/4, for A; e Z. So, the 
continuous-time quantum walk on K2 has the uniform mixing property. We will revisit this analysis under 
more general conditions in a subsequent section. 



4 Circulant Graphs 

A matrix A is circulant if its A;-th row is obtained from the 0-th row by k consecutive right-rotations. A 
graph G is circulant if its adjacency matrix is a circulant matrix. Some examples of circulant graphs include 
complete graphs and full-cycles. An important property of circulant matrices is that they are (unitarily) 
diagonalizable by the Fourier matrix 

F = (14) 



n 



where ui = e^'^*/" and V{ijj) is the Vandermonde matrix given by 



/I 1 

1 UJ 

1 



UJ 



UJ 



UJ 



n-1 
2(n-l) 



(15) 



Let the j-th column vector of V{uj) be denoted by \ujj), j = 0, . . . , n — 1. It is easy to verify that F is 
unitary, i.e., F^^ = F"^, since the Vandermonde matrix obeys V{uj)~^ = V{u>^^). 
If C is a circulant matrix whose 0th column vector is |/) then 



FCF^ = diag{Viu;)\f)). 



(16) 



This Equation El shows that the eigenvalues of a circulant matrix can be obtained by applying the Fourier 
transform F to the first column vector of C. 
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5 Classical versus quantum walks 



5.1 The complete graph Kn 

The adjacency matrix of is A = Jn — In, where J„ is the all-one nxn matrix and I„ is the n x n identity 
matrix. The eigenvalues of i^^A are 1 (once) and — (n — 1 times). Using the orthonormal eigenvectors 

\Fj) = -^\<^j) (the columns of the Fourier matrix F), the initial amplitude vector is |^o) = ^ We 
get 

{— it(n — 2) / \ 

_2g 2(n-i) sinl,J^) ifjVO 
i(e-** + (n - l)e**/("-i)) if j = 



Thus 



1 ^ i ' sm 



Notice that to achieve uniformity, for all j ^ 0, -Pt(O) = Pt{j) or 

4 r, f tn \ 

-sin^ TT, = 1- (19) 

77, \2(n — 1)/ 

Uniformity is possible only if n = 2, 3, 4. In contrast, a classical walk never achieves uniform on K2 but 
converges to uniform for all Kn, n > 2. This is summarized in the following theorem. 



Theorem 1 No complete graph, except for K2, K^, and K4, has the uniform mixing property under the 
continuous-time quantum model. 

5.2 The balanced complete multipartite graphs 

Let G be a complete a-partite graph where each partition has 6 > 1 vertices (the case 6 = 1 is the complete 
graph case). Let A be the normalized adjacency matrix of G. Note that A is given by 

A = -^Ka®\jb. (20) 
a — 1 

Using the circulant structure of both matrices, it is clear that the normalized eigenvalues of Ka are 1 (once) 
and — (with multiplicity a — 1) and that the normalized eigenvalues of Ji, are 1 (once) and (with 

multiplicity 6—1). Let a = e~ be the principal a-th root of unity and let /3 = e~5" be the principal 6-th 
root of unity. Both matrices have the columns of the Fourier matrix of dimensions a and h, respectively, as 
their orthonormal set of eigenvectors. So let { | ctj ) }^^q and {\f3k) }fc=o orthogonal sets of eigenvectors 
of Ka and J^, respectively. Note that |ao) = |la) and |/?o) = jib). 
The eigenvalues of A is Aq = 1 (once) with the all-one eigenvector 

\la)®\h), (21) 

Ai = — (with multiplicity a — 1) with the eigenvectors (with lengths ah) 

j = l,2,...,a- 1. (22) 
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and A2 = (with multiplicity a(6 — 1)) with eigenvectors (with lengths ab) 

® \Pk), I <j <a, 1 <k <b. 
Thus, the wave amplitude function \ipt) is given by 

a— I b—1 a—1 



1 

ab 



e-'*|la) + e-1 \aj) (^\h) + Y,Yl '"^^^ ® '^'^ 

i=l k=l j=l 



which equals to 



1 

a6 



The wave amplitude expressions are 



+ 



\0b) 



' ^(e-^* + e^(a-l) + a(6-l)) ifj = 

J_(^-it 
ab 

y ab\^ 



a 



if 1 < J < 6 
otherwise 



Forcing the uniformity on the last type, for j > b, one obtains 



Pt{j) 



{ab)- 



■ sm 



ta 



2{a - 1) 



1 

a6' 



(23) 



(24) 



(25) 



(26) 



(27) 



which yields the condition sin (ta/[2(a — 1)]) = ^ or 1 < a6 < 4. There are only three legal (integral) 
cases, namely, a = 1, 6 = 4 (or K4), a = 6 = 2 (or K2^2), and a = 4, 6 = 1 (or K/C). The case for the empty 
graph Ki is obvious and we know from Theorem [Qthat has uniform mixing. For the case of i^2,2> the 
three amplitude expressions are |(e** + e~** + 2) = i(cos(t) + 1), |(e** + e~** — 2) = i(cos(t) — 1), and 
I sin(t). For t being odd multiples of 7r/2, their probability forms achieve uniformity. So 7^^2,2 has uniform 
mixing. 

Theorem 2 No balanced complete multipartite graph, except for i^'2,2. ^<^^ the uniform mixing property 
under the continuous-time quantum walk model. 



6 Other graphs 

6.1 The complete cycle C„ 

The adjacency matrix A of Cn is given by 



A 



/ 1 

1 1 

10 

V 1 



1 \ 





1 0/ 



(28) 
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Let H = be the normalized adjacency matrix of the full-cycle C„ on n vertices. Let ujj = uo^ = e^'^*-'/", 
for j = 0, 1, . . . , n — 1. Using the properties of circulants, the eigenvalues of H is given by 

A, = ^(w,- + uo"^-^) = cos(27rj7n). (29) 

The wave equation gives us 

n-1 

n ^-^ n ^-^ 

3 3=0 

Then 

^ n— 1 

(A;|Vt) = - y e-^*™'^(2^-'V")a;f . (31) 

3=0 

This is a complicated nested exponential sum^ Despite the fact that C3 (which is equivalent to 7^3) and C4 
(which is equivalent to K2^2) both have the uniform mixing property, the general case of Cn, n > A, remains 
intractable. 

Conjecture 1 No complete cycle Cn, except for C-^ , C^i, has the instantaneous uniform mixing property 
under the continuous-time quantum walk model. 

Note that C„ has the uniform mixing property in the discrete quantum walk model f^. 



^ n— 1 

1 y^g-itcos(27rj/n) 

n 



3=0 



(30) 



6.2 The Cayley graph Sn 

Let Sn be the symmetric group on n elements (the group of all permutations on [n]) and let T„ be the set 
of transpositions on [n]. The Cayley graph Xn{Sn,Tn) is defined on the vertex set Sn and (7r,r7r) is an 
edge, for all i: € Sn and r G T„; it is a bipartite, connected (2) -regular graph on n\ vertices. Consider 
the simplest case of S-^ with elements {e, (1 2), (1 3), (2 3), (1 2 3), (1 3 2)}. This Cayley graph on ^3 is 
isomorphic to -ftr3 3. By Theorem |2j the continuous-time quantum walk on it is not uniform mixing. We 
have verified that the same is true of 5*4. 

Claim 3 The Cayley graph of 84^ is not instantaneously uniform mixing. 

We conjecture that this phenomenon holds for all n > 4. 

Conjecture 2 No Cayley graph Xn, except for X3, has the instantaneous uniform mixing property under 
the continuous-time quantum walk model. 

It is likely that tools from group representations and characters might settle this question. This is because 
the characters of 5„ yield information about the eigenvalues of the Cayley graph. Given this, one could 
build the set of orthonormal eigenvectors. Unfortunately, there is no known explicit formula that gives all 
the characters of Sn', only estimates are known for the general case. 

'An alternative expression is (fc|V't) = 'n,v=±k (mod N)i~^)'' '^"i^)' where Jy(t) is the Bessel function 1121 . 
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